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Local or nonlocal character of quantum states can be quantified and is subject to various bounds 
that can be formulated as complementarity relations. Here, we investigate the local vs. nonlocal 
character of pure three-qubit states by a four-way interferometer. The complete entanglement in 
the system can be measured as the entanglement of a specific qubit with the subsystem consisting 
of the other two qubits. The quantitative complementarity relations are verified experimentally in 
an NMR quantum information processor. 
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I. INTRODUCTION 

Classical physics groups physical entities in categories 
that are considered to be mutually exclusive, such as 
waves and particles. However, many experimental re- 
sults are incompatible with this approach, since they can 
not be explained in terms of a pure wave picture or a 
pure particle picture. Quantum mechanics was devel- 
oped to resolve this discrepancy and Bohr introduced the 
concept of complementarity [I] to emphasize the different 
approach. The most familiar aspect of complementar- 
ity is perhaps the wave-particle duality. It means, e.g., 
that light has characteristic properties that are usually 
associated with particles but also show behavior usually 
associated with waves. If we design an experiment to 
measure any of these properties, it can only be achieved 
at the cost of losing information about the other. 

Complementarity is often illustrated by means of a 
two-way interferometer such as Young's double-slit ex- 
periment or a Mach-Zehnder setup. Already in 1909, the 
interference pattern of "single photons" was observed ex- 
perimentally in a double-slit interference experiment by 
Taylor [2] and later by Dempster and Batho [3]. This 
was of course possible only because these experiments 
did not provide any information about the path taken by 
the photons in the double-slit interferometer. The same 
effect was also observed with many other kinds of single 
quantum objects including electrons [I] [5], neutrons [6], 
trapped ions [7], atoms [5], and even molecules [5]- 

At the qualitative level, complementarity is thus a well 
established concept. More recently, it was found that 
complementarity can be quantified [10J QTJ Q2J HU HU [15] 
16 . For the case of the wave-particle complementarity, 
it is possible to formulate it in terms of the inequality 

P 2 + V 2 <1. (1) 
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In this expression, the particle-like property is quantified 
by the predictability P, which specifies a priori knowl- 
edge of the path that the system will follow ("which- 
way" information), whereas the wave- like properties are 
quantified by the visibility V of the interference fringes. 
In the case of a pure quantum state, the inequality turns 
into the limiting equality. While this wave-particle dual- 
ity was considered mostly in two-path interferometers, it 
is possible to generalize it to multi-path interferometers 

Quantitative complementarity relations exist not only 
for individual quantum systems, but even more for com- 
posite systems. In systems consisting of two quantons, 
some new complementarity relations were found, such 
as the complementarity relation between single and two- 
particle fringe visibilities [TH [18], between distinguisha- 
bility and visibility[15 , and between the coherence and 
predictability [IB] in a quantum eraser [T9"]. These proper- 
ties are less directly measurable, but some can be quan- 
tified, e.g. by two-particle interferometry. Many of these 
complementarity relations have been experimentally in- 
vestigated by interferometric experiments, using a wide 
range of composite two-quanton systems including pho- 
tons [201 EH [221 E3l El], atoms [25] [26] and nuclear spins 
in a bulk ensemble [27112^11^9] . 

In the course of the study of complementarity in com- 
posite systems, entanglement is found to be a key en- 
try. As a purely quantum correlation with no classi- 
cal counterpart, entanglement can be used to quantify 
the non-local aspects of the composite system. Some 
progress has been achieved in this direction, such as the 
complementarity relations between distinguishability and 
entanglement |3"U] . between spatial coherence of biphoton 
wave functions and entanglement [3T] , between local and 
nonlocal information [32 , and a beautiful equality be- 
tween visibility, predictability and entanglement in pure 
two-qubit states [33 . Additionally, some complementar- 
ity relations in n-qubit pure systems were found, such as 
the relationship between multipartite entanglement and 
mixedness for special classes of n-qubit systems [34], and 
between the single particle properties and the n bipartite 
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entanglements in an arbitrary pure state of n qubits [3"5]. 

In our previous paper |28j . we found a complementarity 
relation that exists in an n-qubit pure state: 



II. EXPERIMENTAL SETUP FOR A 
THREE-QUBIT SYSTEM 

A. Preferred basis 



c 



k(ij...m) + ^fe 



si = 1. 



(2) 



This relation implies a tradeoff between the local single- 
particle property (S 2 ) whose two constituents are P% and 
V fe 2 , and the nonlocal bipartite entanglement between the 
particle and the remainder of the system (C^y m ^ ) , de- 
fined in terms of the marginal density operator p k [57] 



C, 



k{ij. ..m) 



2[1 



Tr(pl)}. 



(3) 



Moreover, a conjecture was made: the bipartite entan- 
glement C?, - ■ , might be equal to the sum of all possi- 
ble pure multi-particle entanglement (s) connected to this 
particle [25J. This conjecture was proved for pure two- 
and three-qubit systems |28j . Therefore, measuring the 
bipartite entanglement C 2 ^- m ^ implies that we obtain 
an entire entanglement (nonlocal) connected to this par- 
ticle. Therein the simplest case with two qubits has been 
verified by NMR interferometry, i.e., C 2 + S 2 = 1, where 
C is the concurrence of a two-qubit state ip which is re- 
lated to "the entanglement of formation" [35], defined 

by 



CO) = \ {lp\(J V ® 0-y\ll)* 



where a y is the y component of the Pauli operator and 
\ip*) is the complex conjugate of l^). 

The question that was left open in this earlier paper 
is, if it is possible to test the complementarity relation 
in a system with more than two qubits. The present 
paper shows an example for such an experimental test 
in a pure three-qubit system. For a pure state |£) of a 
three-qubit system ABC, we use a generalized four- way 
interferometer to verify the complementarity relation 



c 2 A(BC) (\0) + s 2 A (\0) = i- 



(4) 



This experiment uses a specific property of pure three- 
qubit states. In the next section, we will describe this 
property and the experimental configuration used to 
measure the quantities of Eq Q. Sec. Ill adds details 
about the main components (transducers) in the inter- 
ference experiment. In section IV, we combine the inter- 
ferometer with state preparation and readout. Section 
V gives experimental details of the implementation in an 
NMR quantum information processor for different classes 
of pure 3-qubit states and discusses the results. 



Let us express the pure state |£) of the three-qubit 
system ABC in the standard basis \ijk) k — 0, 1): 



|0 = J2a ljk \ijk). 



(5) 



The coefficients aijk are normalized to 1. If we regard the 
pair BC as a single object, it makes sense to consider the 
concurrence Ca(bc) between qubit A and the composite 
object consisting of the two qubits B and C. 

An interesting and unique property of a pure state |£) 
of the three-qubit system helps us to design an experi- 
mental scheme for measuring the concurrence Ca(bc) by 
an interference experiment: The reduced density matrix 
Pbc has at most two nonzero eigenvalues. Accordingly, 
even though the state space of BC is four dimensional, 
only two of those dimensions are necessary to express the 
state |£) of ABC [35]. Therefore, the state |£) can always 
be rewritten as 



io = em^: 



BC 



(6) 



where are the eigenstates with the two nonzero eigen- 
values of the reduced density matrix psc = TYa(|£) (0); 
and the real coefficients bij are normalized to 1. There- 
fore, we can treat A and BC, at least for the present 
purpose, as a pair of qubits in a pure state. 

Like in a two-qubit system, we can thus design a four- 
way interferometer for a pure 3-qubit state, which con- 
sists of four paths, which we label by the correspond- 
ing basis states {\0)a, |$o)bc, |$i)bc}- Figure [l] 
shows the reference setup: The source S emits three 
particles A, B and C in a pure state. Particle A can 
propagate along path |0) ^4 and/or through a vari- 

able phase shifter ip\. Beamsplitter BSi connects the 
two paths and the particles are then registered in either 
beam l-f^i)^ or \L{)a- One the other side the pair of 
particles B and C as a whole can propagate along the 
paths \&o)bc and/or \&i)bc- The probabilities of the 
joint (e.g. p(\K 1 ) a \K 2 )bc)) and single (e.g. p(\Ki)a)) 
events generate interference patterns as a function of the 
phase angles ipi, depending on the state |£) of the source 
S. 

From the resulting one-party interference pattern of 
qubit A, the single-particle fringe visibility is defined as 



(7) 



where x = K\ or L\, and p m i n and p max are the mini- 
mal and maximal probabilities (as a function of ipi). The 
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FIG. 1: Schematic four- way interferometer in a pure 3-qubit 
system, using beam splitters BSi, BS2 and phase shifters tpi, 

V?2- 



other quantity related to the single-party property of par- 
ticle A is the predictability Va, which quantifies the a 
priori which-way knowledge. It is defined as 



(8) 



where a, = 



1 

-1 



is the z component of the Pauli 



oprator. Va thus measures the magnitude of the prob- 
ability difference that particle A takes path |0) or the 
other path 

We combine these two single-party properties into a 
single entity 



v\. 



which measures the single-particle character for particle 
A. 

The two-party (nonlocal) properties between qubit A 
and the pair of qubits B and C can be measured by higher 
order correlations. Following references [TH[T5], we use 
the "corrected" two-party fringe visibility 

V [PdOlmax- [P(IO)] min , Q x 

Vw(l0) ^ii^m:' (9) 

where the state |£) is the product state \x) A \u)bc with 
x = K\ or L\ and y = K% or L^. The "corrected" joint 
probability p is defined as 



p{\x) a \v)bc) =p{\x) a \v)bc) ~ P (\x) A ) P (\v) bc) + 



1 



This correction eliminates single-party contributions [14 

KB]. 

The single-party and two-party properties satisfy a du- 
ality relation: 



V, 



A(BC) 



si = 1. 



(10) 



Here |ifi) A , \Li)a is an arbitrary basis in the Hilbcrt 
space Ha of particle A. To get the equality, we have to 
choose a specific basis for the BC subsystem: \K2)bc, 
\L2)bc must be linear combinations of the two states 
that correspond to the nonzero eigenvalues of the reduced 
density operator pbc- We will refer to this basis as the 



preferred basis. With this basis, the two-party visibility 
Va(bc) becomes equal to the concurrence Ca(bc)i i.e., 



v. 



A{BC) 



c 



A{BC) 



2|&oo&ii - b Q1 b 10 \ 



This was proved in our previous paper [55] • Therefore, 
the concurrence Ca(bc) °f the source S can be quantita- 
tively measured by the two-party fringe visibility Va(bc), 
so as to verify the complementarity relation (1). 



B. Extended basis 

For most of the calculation we assume that the 
measurement basis for the BC subsystem consists of 
two states that are within the subspace spanned by 
{|$o) bc, |$i) bc}- It is possible to choose a different ba- 
sis, and, for an unknown input state, it is not possible to 
choose a basis that falls into the {\&o)bc, \®i)bc} sub- 
space. In the general case, the paths in the BC part of 



the interferometer must be written as 



(with normalized coefficients Cj). 
probability is then 



The "corrected" joint 



p(\0) a Hbc) 



i=0 



|2= 



P(\0) A \^ 



•BC) 



and the two-party visibility becomes 
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Va(bc) (\0) a Hbc) = £ \ci\ 2 V$ BC) 



with V 



(i) 



A(BC) 



and < V 



(i) 



A(BC) 

{V$b C) } < V 



Va(bc) (|0) a \*i) B c)- Since < \c t \ 2 < 1 
< 1, we find 



A(BC) 



(10), 



\™)bc)< max { V l 4 (BC)J'' 



V 



A(BC) 



si <i. 



(ii) 



The limiting case of the limiting equality ( 10 1 is obtained 
if two conditions are fulfilled: (i) the preferred basis is 
chosen as the measurement basis and (ii) the transducer 
acting on the BC subsystem acts only on the subspace 
{|S )bc,|S 1 )bc}- 



III. TRANSDUCERS 

A. Preferred basis 

Both parts of our interferometer (Fig. 1) contain a 
transducer consisting of a variable phase and a symmet- 
ric beam splitter. As discussed in Ref. [25], this combi- 
nation provides a universal interferometer. Mathemati- 
cally, they can be described by the unitary operation hi, 
written in the preferred basis {|j)yi| < i ) i)sc} : 



U (<pi,<p2) = Ua (ipi) <8> Ubc (<pz) 



(12) 
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Each transducer Ua ('Pi) and Use (^2) maps the input 
state into an output state by the transformation: 

1 / e -i<Pi/2 e i Vi /2 \ 

U ^ = V -e-^'/ 2 e^/ 2 ) ■ (13) 

Here we use Ubc {^2) to represent the matrix expression 
in the preferred basis {|$o) bc> l^i) bc}- 

B. Extended basis 

However, although we can regard the pair BC as a fic- 
titious qubit spanned by the vectors {\&o)bc, \®i)bc}, 



the practical operation in the experiments on the ob- 
ject BC is four dimensional. This requires us to con- 
struct a four-dimensional unitary operation Ubc m an 
orthonormal basis {I'&o), l^), ^3)}, m whose sub- 
space {|$o);|^i)} the transformation has the effect of 
Ubc (^2)1 while it acts as an arbitrary single qubit oper- 
ator on the subspace {l^), l^a)}- Therefore, the trans- 
formation Ubc in the basis {|<&o), ^2), ($3)} can be 
written in the form 



U B c{V2) = 



where a, (3,^,5 are real numbers. A relative simple way is 



1 

V2 
1 



V2/2 J_ p ^V2/2 

i V2 /2 J_ 
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71 
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— e -i<P2/2 
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with 7= § and (5= -8 = 



e iv>2/2 
e i<P2/2 






-iS 



- sin 

cos |e I/3 



(14) 






e -iip 2 /2 
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J<p 2 /2 
0W2/2 



(15) 



IV. NETWORK 
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FIG. 2: Schematic network for the experimental verification 



of the complementarity relation C A ^ BC ^ + S A = 1 in a pure 
3-qubit system. The pseudo-Hardmard gate H — e l 4 ay ro- 
tating the qubit by the angle J about the —y axis, the phase 

• ^1 

shift gate Z(ipi) — e~'^ <T ^ rotating the qubit by the angle tpi 
about the z axis, and the gates 1Z and Ubc (^2) are explained 
in the text. The last two operations are inverses of each other 
and can be omitted in the experiment. 

Fig. [2] shows the network corresponding to the inter- 
ferometer of Fig. [T] for a pure 3-qubit state. After the 
preparation of a pure 3-particle source, the transducer 



Ubc (^2) must be realized. To do this, we can use the 
notation: 

Ubc{V2) = K>~ 1 Ubc (y>a) ft 

where the operation 1Z transforms the chosen basis of the 
interferometer (e.g., here {\®o}, |$x)> I ^2)1 1*^)}) to the 
computational basis {|00), |01), jlO), |11)}. 

The measurement observable is defined in the pre- 
ferred basis {|$o),|^i)}, while the experimental de- 
tection scheme operates in the computational basis 
{ 1 00) , 1 1 ) , 1 1 0) , 1 1 1 > } . The dashed box labeled "Measure- 
ment" in Fig. [2] therefore starts with a basis transforma- 
tion 1Z, which is followed by the projective measurement 
in the computational basis. 



V. EXPERIMENTAL TEST 
A. System 

As a quantum register for these experiments, we se- 
lected the three 19 F nuclear spins of Iodotrifluoroethylene 
(F 2 C=CFI), shown in the inset of Fig. [3] This system 
has relatively strong couplings between the nuclear spins, 
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large chemical shifts, and long decoherence times. The If we choose a 
Hamiltonian of this system is (in angular frequency units) a W _ q/ 11 ) 



(0) 



„(oo) 



,,(01) 



,,(10) 



7r, we obtain a GHZ-class state 



and 



(16) 



i<3 



where the /* are the local spin operators. Qubits 1, 2 
and 3 represent particles A, B and C of section |Tlj 

Fig. [3 shows the 19 F NMR spectrum of this molecule, 
together with the relevant coupling constants. The lower 
part contains the full spectrum, with the groups of lines 
labeled by the index of the qubits. The upper part shows 
the partial spectra of each qubit on an expanded scale. 
Each qubit is coupled to the other 2 qubits, resulting 
in four resonance lines. In the figure, we have labeled 
these lines with the corresponding logical states of the 
coupled qubits. The numerical values of the coupling 
constants Jjj are given in the inset, together with the 
molecular structure. The relaxation times are T\ = 5.6 s 
and T 2 = 1.9 s. 

The experiments were performed on a Bruker Avance 
II 500 MHz (11.7 Tesla) spectrometer equipped with a 
QXI probe with a pulsed field gradient. The resonance 
frequency for the 19 F spins is « 470.69 MHz. 



B. Initialization 

In the NMR experiments, the system was first pre- 
pared in a pseudopure state (PPS) [301 E] Pooo = 
^jpl + e 1 000) (000 1 with 1 representing the unity oper- 
ator and e w 10~ 5 the polarization, instead of the pure 
state 1 000). Starting from thermal equilibrium, we used 
spatial averaging [42] to prepare the PPS; the pulse se- 
quence [43] is shown in the first part of Fig. [4] 

In the whole experiment, we used robust strongly mod- 
ulating pulses [331 [351 US] to implement all local gates 
(e.g., [f]_ , and [tt]*' 2 etc.). In order to confirm the 
state preparation, we performed a complete state to- 
mography |47j to reconstruct the experimentally normal- 
ized relevant pure part p exp of the density matrix p: 
p exp = p — ^jpl, which is shown in Fig. pi The ex- 
perimentally determined state fidelity [3E] was 



F(p t h,Pexp) = Tr{^J s/pthPexp-s 



0.99. 



With respect to scale-indenpent NMR observations and 
unitary evolution, a pseudo-pure state is equivalent to 
the corresponding pure state [I9"l 150 1 I5T 1 |52] |55j |551 1551 
I5B1 15?] . Therefore, we focus on the relevant pure part 
Pexp of the pseudo-pure state in the remaining sections. 



Oil 



lOcffz = cos -HOOO) 



■sinylm). 



(17) 



The corresponding pulse sequence is represented in the 
second part of Fig. [4j This type of states contains only 
tripartite entanglement, no bipartite entanglement. 



.,(!) 



(01) 



With a' 2 ' = ctg ' = a 3 
we obtain a W-class state 



(10) 



a 



(11) 



and a 



(00) 



10 



w 



(0) 

cosy cos^-|001) 

(0) 

cosy sin -|- 1 010) - 



sin ^|100). (18) 



These states contain only bipartite entanglement, but no 
tripartite entanglement. 



For a. 



(1) _ (01) _ (10) _ (11) _ 



0, we obtain the 



state 



\0i 



„(oo) 



cos — COS — — COS — ^ 
2 2 2 



-|000) 



Oil 



.,(0) 



,,(00) 



cos — cos sm 

2 2 2 

(0) 

■ cos y sin ^- 1 010) 



|001) 
sin ^|100). (19) 



It has tripartite as well as bipartite entanglement. These 
states are intermediate between the GHZ-class and the 
W-class states. We call them the intermediate-class 
states. 



D. Implementation of the Interferometer 

To implement the interferometer shown in Fig. [T] 
we first need to determine the relevant eigenbasis 
{I'&o), | < 1 ) i)}bc of the reduced density operator psc of 
the BC subsystem. Table [I] lists the eigenbases for the 
three initial states that we consider in this context. 

To apply the transducer in this basis, we have to find 
the basis transformation between this eigenbasis and the 
computational basis. Writing |$q) as a general two-qubit 
state, 



C. Preparation of the 3-qubit source 

The 3-qubit states |£) = J2i j k a ijk\ijk) of Eq. (j5j) 
were prepared from |000) using the quantum circuit of 
Fig. [6j The resulting coefficients are 



cos( 



Oil 



(*) 

2*)cos( — 



^j)cos(^ -k). 

2 2 2 



|*o) = cosycosy|00) + sin^sin^|01) 

+ sin ^ cos 1 1 10) + cos ^ sin \ 11) , (20) 

we can transform it into the computational basis state 
1 00) by the transformation 
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FIG. 3: 19 F NMR spectrum of the Iodotrifluoroethylene, measured in a field of 11.7 T. The chemical shifts are 12020 Hz, Hz 
and -17330Hz. The inset shows the structure of the molecule and the relevant coupling constants. 



TABLE I: The preferred basis {|<&o}, |$i)}sc in the interferometer (see Fig. [T]) for three different classes of input states. For the 
intermediate-class states, the coefficients are a\ = cos = sin sin = sin cos ^ for |$o) an d a 'i = — sm ^r, 012 = 

cos sin e -f, a' 3 = cos cos ^ for |$ x ). 



State class 


l*o) 


1*1) 


GHZ 
W 

Intermediate 


|00) 
|00) 

ai| 00) +o 2 |01) +o 3 |10) 


|H) 

(0) (0) 

cos ^-| 01) + sin^§-|10) 
oi|00)+o' 2 |01} + 4|10) 



/ cos ^ cos ^ sin ^ sin ^ sin \ cos ^ cos \ sin ^ \ 



n = 



cos y sm 2 si n 2 cos 2 
cos ^ sin & 



— sin -J- cos 



sin 2 sin f 



cos 



COS 



- cos ^ sin ^ 



COS -y COS y 



(21) 



The same operator also maps the other basis states 
= 1, 2, 3 into basis states of the computational ba- 
sis. For the GHZ states, the pulse sequence for the im- 
plementation of this transformation is shown in the third 
part of Fig. [IJ Depending on the basis states {|$»)}, a 
permutation of the computational basis states is required. 

The preferred bases for the different states are shown 
in Table [T] For the GHZ-class states, the parameters are 
9i = 0, (i = 1,2,3) and the transformation operator 1Z 
simplifies to CNOT 32 . For the W-class states, the param- 
eters are 61 = 62 = 0,63 = — al . For the intermediate- 



class states, the parameters Qi are determined by the the 
two non-zero eigenvalues of the reduced density matrix 
Pbc, 



where 



X± = \(l± Jl-4Asin 2 ^, 



„(0) „(00) 



.,(0) 



A = cos 2 ^(cos 2 ^- sin 2 ^— + sin 2 "M. 

2 y 2 2 2 ' 

In terms of these parameters, we find for the parameters 
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FIG. 4: (Color online) Sequence of radio-frequency and field gradient pulses (PFG) used to prepare the pseudopure initial state 
1 000} (first part), transform it into a GHZ-class state (second part), implement the interferometer (third part), and measure 
the resulting interference patterns (fourth part). The phases of the n pulses in the interferometer, applied to qubits 1 and 2, 
are rji = 7r ~' Pi . [a]^ denotes a pure bilinear evolution of qubits i and j for a duration 2jr " , which is realized by the refocusing 
scheme shown at the bottom. 
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FIG. 5: (Color online) Experimentally measured relevant pure 
part p e:L p of the prepared pseudo-pure state pooo reconstructed 
by tomography. The rows and columns are numbered with the 
computational basis states. 
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FIG. 6: Quantum circuit to prepare a pure 3-particle source 
state \Qabc- Horizontal lines represent qubits; time runs 
from left to right. Conditionality on the other qubits being 
in the |1) and |0) state are represented by filled and empty 
circles, respectively. 



and 



a AO) (00) 

"3 C-i . a i 

tan — = cot — — sin — ^ — . 
2 2 2 



E. Measurement 



9i in the basis transformation 1Z 



tan 2 -i 



X.-A' 



After the interferometer, the output state of the com- 
plete 3-qubit system is 



Wo 



Ua^SK-WUbcMPKIO, (22) 



where P is the relevant permutation matrix. For these 
states, we measure the joint probabilities p{\i) j) bc) 



for detecting particle A on port |z) and particles B and C 
on port of the interferometer. This probability can 
be written as the projection of the output states onto the 
measurement basis 



p(\i) A \$j)Bc) = |sc($jU(j|Vw)| 5 



(23) 



This expression can be evaluated in the computational 
basis \ijk) by using the transformation operator 1Z: 

p{\i) A \*j)Bc) = \B C m A {i\R\^ out )\ 2 . (24) 
The single-particle probabilities of particle A are 

p{\i) A )= A (A Tr Bc(\lpout)(lpout\)\i}A- (25) 

|10>I2|00» 12 |H»12|01>12 |11>13 |M>,3 |11>, 3 |01> 13 |11> 2 3 |l»>23l»l>23 l»»>23 
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FIG. 7: Experimental spectra of the three 19 F spins (corre- 
sponding to qubits 1, 2 and 3 from right to left) for the GHZ 
state when tp = tpi — tp2 = —n/2. 

To determine the probabilities, we measured the pop- 
ulations of all eight computational basis states by first 
deleting coherences with a field gradient pulse and then 
applying selective readout pulses to the individual qubits. 
This procedure is denoted by the last part of Fig. [4] The 
dashed read-out pulses indicate that the three pulses were 
applied in three separate experiments. From each of the 
three FIDs, we obtain a spectrum of the corresponding 
spin with four resonance lines after Fourier transforma- 
tion. Fig. [7] show representative spectra for the case of 
the GHZ state and a phase shift of tp — p\ — p2 = —tt/2 
in the interferometer. 

As the path length of the interferometer arms is 
changed, interference between the two paths changes the 
populations of the different states. This oscillation can 
be observed directly in the NMR spectra, as shown in 
Fig. [5J where we have plotted the variation of the sub- 
spectrum of qubit 1 as a function of the phase ip for the 
GHZ state. 

The relevant populations were obtained from the spec- 
tra by integrating over the resonance lines. Fig. [9] shows 
interferograms for some of the populations (equal to the 
probabilities p{\i) A \^j) bc) for the GHZ-state (left hand 
side) and a product state (right hand side). Clearly, the 
maximally entangled state shows high visibility fringes, 
while the variation of the populations essentially vanishes 
for the product state. 

We extracted the visibility V A (bc) from the interfer- 
ograms of the populations p(\i) A \<& j) bc) by Eq. (pi. 
The single-particle probabilities p(\i) A ) were obtained by 
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FIG. 8: Variation of the spectrum of qubit 1 with the inter- 
ferometer phase tp for the GHZ state. 



summing the the populations related to the state \i) A , 
from which we obtained the singe-particle visibility V A . 
The results is summarized in the next subsection. 



GHZ state 



Product state 



0.5 



a 0.25 
o 
P. 







o 


p(l0>l* >)=p(l000>) 






+ 


p(io>id>,>)=p(ion>) 






V 


p(ll>l<J>„>)=p(H00>) 






4- 


plll^^pCIlll^ 




n So n 























































-it/2 

Phase angle qp 



jc/2 -jt/2 re/2 

Phase angle rp 



FIG. 9: (Color online) Relevant joint populations recon- 
structed from the experimental data for (a) the GHZ state 
and (b) the product state. 

In addition to the visibility, we measured the pre- 
dictability V A . According to Eq. Q, it is given as the 

expectation value of the observable ai A \ i.e. the popula- 
tion difference of particle A. We measured it by applying 
a field gradient pulse that destroyed the coherences and a 
subsequent readout pulse [^] that converted into 

ai A \ which was then recorded as the FID. After Fourier 
transformation, we integrated over the relevant resonance 
lines in the spectrum. The absolute value of this integral 
corresponds to the predictability V A . 



F. Experimental results 



Fig. 10 summarizes the experimental results for three 
different classes of states by plotting the degree of local 
vs. nonlocal character for each case. According to sec- 
tion II, they should be related by the complementarity 
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relation V^bc) 
where Sa 



1 for the pure three-qubit states, 



Pi 



Va(bc) the nonlocal character 



A quantifies the local character and 
Clearly, the experimen- 
tal data points (circles) agree well with the theoretical 
prediction (solid curves). In the first system (GHZ), the 
non-local character is due exclusively to tripartite entan- 
glement; in the second system (W), it arises from bi- 
partite entanglement, and in the third case, we have a 
combination of both types of entanglement. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

Nonlocal character \a(BC) 



FIG. 10: (Color online) Experimental verification of the com- 
plementarity relation of V\< B c\ + <Sa = 1 in a pure three-qubit 



system: (a) for GHZ-class states of Eq. (171, (b) for W-class 
state s of Eq. ( 18 1 and (c) for intermediate-class states of 
Eq. (19 1. Solid curves represent the theoretical complemen- 



tarity relation of local character Sa versus nonlocal property 
Va(bc) ■ Experimental results are indicated by circles. 

The deviation between the experimental and theoret- 
ical values is primarily due to the inhomogeneity of the 
radio frequency field and the static magnetic field, imper- 
fect calibration of radio frequency pulses, and signal de- 
cay during the experiments. The experimental errors are 
bigger for states with predominantly "non-local" charac- 
ter. This is well compatible with the general expectation 
that nonlocal states are less robust and is observed even 
for the initial state preparation, where the measured fi- 
delity is lower for the nonlocal states. For example, the 
experimental fidelity of the GHZ state ^(|000) + 1 1 1 1> ) 



with ai = tt/2 in Eq. (17 1 is about 0.97, compared to 
the 0.99 fidelity of the product state |000) with a x = 
inEq. (pi 



formulated. Here, we have quantitatively compared the 
local versus nonlocal character of the quantum states of 
three coupled qubits (spins 1/2). 

For the experimental measurements, we used a sys- 
tem of three coupled nuclear spins in a liquid-state NMR 
spectrometer. The degree of local vs. nonlocal char- 
acter was measured by constructing a suitable four-way 
interferometer and utilizing a specific property of pure 
three-qubit states: In any two-qubit subsystem psci a t 
most two eigenvectors of the density matrix have non- 
zero eigenvalues. This allowed us to quantify the local 
character of the system by measuring the polarization of 
the (arbitrary) particle A and the nonlocal character via 
a measurement of the entanglement between the single 
qubit A and the subsystem BC . While the interferom- 
eter only uses four channels, they were chosen in such 
a way that the measurement results quantify the com- 
plete entanglement of the three-qubit system, including 
bipartite as well as tripartite contributions. 

Here we have restricted our theoretical treatment to 
cases where the three-qubit system is described in a pure 
quantum state. In an experiment, the prepared states 
inevitably involve some mixture less or more. Naturally, 
the question arises as to what happens when the system 
is initially a mixed state. As discussed in Ref. [33] and 
[2"8] for a bipartite system, a weaker statement for the 
complementarity of Eq. Q is believed in the form of an 
inequality C A( ^ BC ^(p) + S^(p) < 1 for the most general 
case of a mixed three-qubit system. However, similar 
to the case of a bipartite system [251 [33], there is no 
corresponding inequality for the visibility Va(bc) m the 
mixed three-qubit states because the definition of Va(bc) 
by Jeager et al. [TU [TS] is unfeasible and the direct rela- 
tion between the concurrence Ca(bc) an d the visibility 
Va(bc) ceases to exist for mixed states. 

The complementarity relation that we have verified 
here, can be used to measure the degree of entanglement 
by measuring only the single-particle character of a given 
state. Our measurements extend earlier tests of comple- 
mentarity that were done in one- or two-qubit systems 
[21 HI [3 El [71 El [201 E2l [231 E3 [261 [271 EH [291 EH]. This 
experiment may thus be considered as a first step towards 
establishing and testing quantitative complementarity re- 
lations in multi-qubit systems. 



VI. CONCLUSION 
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